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Abstract 

Starting from a (new physics independent) tree level determination of p and 77, 
we perform the Unitarity Triangle analysis in general extensions of the Standard 
Model with arbitrary new physics contributions to loop-mediated processes. Us- 
ing a simple parameterization, we determine the allowed ranges of non-standard 
contributions to |A.F| = 2 processes. Remarkably, the recent measurements from 
B factories allow us to determine with good precision the shape of the Unitarity 
Triangle even in the presence of new physics, and to derive stringent constraints 
on non-standard contributions to |A.F| = 2 processes. Since the present experi- 
mental constraints favour models with Minimal Flavour Violation, we present the 
determination of the Universal Unitarity Triangle that can be defined in this class 
of extensions of the Standard Model. Finally, we perform a combined fit of the 
Unitarity Triangle and of new physics contributions in Minimal Flavour Violation, 
reaching a sensitivity to a new physics scale of about 5 TeV. We also extrapolate all 
these analyses into a "year 2010" scenario for experimental and theoretical inputs 
in the flavour sector. All the results presented in this paper are also available at the 
URL http://www.utfit.org, where they are continuously updated. 



1 Introduction 



With the increasing precision of the experimental results, the Unitarity Triangle (UT) 
analysis shows the impressive success of the CKM picture in describing CP violation 
in the Standard Model (SM). UT parameters have been consistently determined using 
both CP-conserving (IK^/Vy, Am^ and Ama/ Am s ) and CP-violating [ek and sin2/3) 
processes [1]. Additional measurements of several combinations of angles of the UT, 
especially 7 from B — > DK and a from charmless B decays, confirm this picture [1]. 

This success becomes a puzzle once, as a possible solution to the gauge hierarchy 
problem, the SM is considered as an effective theory valid up to energies not much higher 
than the electroweak scale. Indeed, even in the favourable case in which the theory above 
the cutoff is weakly coupled, such as the Minimal Supersymmetric Standard Model, large 
contributions to Flavour Changing Neutral Current (FCNC) and CP-violating processes 
are expected to arise [2], clashing with the large amount of accurate experimental data 
now available on these transitions. This is due to the presence of additional sources 
of flavour and CP violation beyond the CKM matrix (see for example Ref. [3] for the 
supersymmetry case). 

In general, the flavour puzzle admits two classes of possible solutions. The first one 
contains models in which a flavour symmetry is invoked to explain the hierarchy of quark 
masses and mixing angles. These models are based on different theoretical approaches 
(supersymmetry, grand unification, extra dimensions, . . . ) leading to different levels of 
agreement with the data and to different low-energy signals. However, low-energy pro- 
cesses generally receive sizable additional contributions which jeopardize the validity of 
the SM UT analysis (see Ref. [4] for a supersymmetric example). A generalized UT fit 
allowing for the presence of arbitrary New Physics (NP) contributions is therefore very 
useful for model building, since it provides at the same time the allowed ranges for the 
SM CKM parameters and for the NP contributions to |AF| = 2 processes. This will be 
the subject of the first part of the present work (Sees. I2EJ). 

The second class of solutions to the flavour puzzle contains models with Minimal 
Flavour Violation (MFV). The basic idea of MFV is that the only source of flavour vi- 
olation is in the SM Yukawa couplings, so that all FCNC and CP-violating phenomena 
can be expressed in terms of the CKM matrix and the top quark Yukawa coupling [5-7]. 
This leads to strong correlations between different observables, and allows for a detailed 
study of low-energy phenomena. While there are several implementations of MFV in dif- 
ferent contexts (two-Higgs doublet models, supersymmetry [8], extra dimensions [9], . . . ), 
it is possible to perform two very general analyses under the MFV hypothesis. The first 
is the determination of the so-called Universal Unitarity Triangle (UUT) [5], which is a 
UT fit performed using only quantities that are independent of NP contributions within 
MFV models. The second is a simultaneous fit of the UT and of NP contributions in the 
\AF\ = 2 sector. These two analyses will be presented in the second part of this work 
(Sec. EJ), and they serve as the starting point for the study of rare decays and CP violation 
in MFV models [10]. 

Finally, in Sec. El we present the possible future improvements in the above analyses 
by considering a "year 2010" scenario for experimental data and theoretical inputs in the 
flavour sector. 

While to our knowledge the determination of the UUT is presented in this work for 
the first time, several attempts have been previously made in the study of the UT in the 
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presence of NP. Considering only model- independent analyses, in Ref. [11] the case of NP 
contributions to \AB\ = 2 or \AS\ = 2 transitions was analyzed: this corresponds to 
the discussion in Section 0] of the present work. A first version of the present analysis, 
with some experimental constraints missing, was presented in Ref. [12]. Constraints on 
NP in the \AB\ = 2 sector using B physics only were considered in Refs. [13-15]. The 
determination of the UT from tree- level processes only was presented in Ref. [1]. A 
general analysis was recently performed in Ref. [16], but not all available constraints were 
used. With respect to these previous studies, we improve several theoretical aspects and 
perform a simultaneous determination of UT and NP parameters using all the available 
constraints in all sectors. 

A compilation of the experimental and theoretical inputs to our analyses is presented 
in Tabled 



Parameter 


Value 


Gaussian (a) Uniform 

(half-width) 


A 


0.2258 


0.0014 


\V cb \(excl) 


41.3 x 10 -3 


1.0 x 10 -3 1.8 x 10 -3 


\V cb \(mc\.) 


41.6 x 10 -3 


0.7 xlO -3 


|Kb|(excl.) 


38.0 x 10 -4 


2.7 xlO -4 4.7 xlO -4 


|Kfc|(incl.) 


43.9 x 10 -4 


2.0 x 10 -4 2.7 x 10 -4 


Am d [ps -1 ] 


0.501 


0.005 


Am s [ps -1 ] 


> 14.5 at 95% C.L. 


sensitivity 18.3 


F Bs \/b Bs [MeV] 


276 


38 


<• Fb s yj Bb s 


1.24 


0.04 0.06 


B K 


0.79 


0.04 0.09 


£k 


2.280 x 10 -3 


0.013 x 10 -3 - 


Jk [MeV] 


160 


fixed 


Am K [ps -1 ] 


0.5301 xlO" 2 


fixed 


sin 2/3 


0.687 


0.032 


m t [GeV] 


163.8 


3.2 


m b [GeV] 


4.21 


0.08 


m c [GeV] 


1.3 


0.1 


a s {M z ) 


0.119 


0.003 


G F [GeV" 2 ] 


1.16639 xlO -5 


fixed 


raw [GeV] 


80.425 


fixed 


iti b o [GeV] 


5.279 


fixed 


m B n [GeV] 


5.375 


fixed 


m° K [GeV] 


0.4977 


fixed 



Table V. Values of the relevant quantities used in the UT fit. The Gaussian and the flat 
contributions to the uncertainty are given in the third and fourth columns respectively 
(for details on the statistical treatment see Ref. [17]). Several branching ratios and CP 
asymmetries have been used. Their values and errors can be found in Ref. [18] and have 
been updated to Summer 2005. 
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Figure 1: JTie selected region on p-fj plane obtained from the determination of \V u b/V c b\ 
and 7 (using DK final states). Selected regions corresponding to 68% and 95% probability 
are shown, together with 95% probability regions forj and \V U b/V c b\- 



2 Unitarity clock, unitarity hands: 

A model-independent determination of the UT 

Let us first of all discuss the shape of the UT in the presence of arbitrary NP contributions. 
All the available experimental data exclude the possibility of sizable contributions to tree- 
level SM processes, so that extensions of the SM in which NP enters low-energy processes 
at the tree level are strongly disfavoured. We can therefore safely assume in this work 
that NP enters observables in the flavour sector only at the loop level. It is then possible 
to determine two regions in the p — f) plane independently of NP contributions, using only 
tree-level B decays. The CKM elements V u b and V c b are determined using semileptonic 
inclusive and exclusive B decays. The angle 7 is obtained by measuring the phase of V u b 
appearing in the interference between b — > c and b — ■> u transitions to DK final states. 1 

1 We neglect possible NP contributions to D°-D° mixing, since their contribution is expected to be 
well below the present experimental accuracy [19]. In the future, it might become necessary to take them 
into account following Ref. [20] . 



4 



As shown in Fig. ^ it is now lunchtime (~13:35) on Andrzej's unitarity clock [21]. 2 

The results of this analysis, reported in Tab. El can be used as a reference for model- 
building and phenomenology in any extension of the SM with loop-mediated contributions 
to FCNC processes. The present precision is expected to improve considerably in the near 
future, as discussed in Sec. |H1 



UT fit 


- using only \V ub /V cb \ and 7 




SM Solution 


2 nd Solution 


P 


0.18 ± 0.12 


-0.18 ± 0.12 


fj 


0.41 ± 0.05 


-0.41 ± 0.05 


sin 2(3 


0.782 ± 0.065 


-0.641 ± 0.087 


7[°] 


65 ± 18 


-115 ± 18 


«[°] 


87± 15 


-46 ± 15 


2/3 + 7 [°] 


122 ± 13 


-152 ± 13 



Table 2: Results for several UT parameters, obtained using the constraints from \V u b/V c b\ 
and 7 (using DK final states). 

Beyond the Standard Model, one can include the information from other constraints, 
taking into account the effect of NP in a general way. In particular, one has to consider 
two effects: 

• The contribution of new operators in the |AF| = 2 Hamiltonian, which affects 
mixing processes and, as a consequence, the determination of Am rf s , Ek and of the 
angles (3 and a. 

• The effect of NP in the |AF| = 1 Hamiltonian, for all those processes occurring 
through penguin transitions. In our case, this concerns the determination of a from 
charmless B decays and the CP asymmetry in semileptonic B decays A^. 

3 Model-independent constraints on New Physics in 
|AF|=2 transitions 

Our goal in this Section is to use the available experimental information on loop-mediated 
processes to constrain the NP contributions to |AF|=2 transitions. In general, NP models 
introduce a large number of new parameters: flavour changing couplings, short distance 
coefficients and matrix elements of new local operators. The specific list and the actual 
values of these parameters can only be determined within a given model. Nevertheless, 
each of the mixing processes listed in Tab.|3]is described by a single amplitude and can be 
parameterized, without loss of generality, in terms of two parameters, which quantify the 
difference of the complex amplitude with respect to the SM one [22]. Thus, for instance, 

2 Fig- HI first appeared in Ref. [1]. Similar results were recently obtained in Ref. [16]. 
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Tree-Level B° mixing K° mixing B° mixing 

\V ub /V cb \ Am d e K Am s 

7 (DK) A CP (B - J/^K) A CP (B° S - J/^0) 

A C p(B -> nn,pir,pp) 



Table 3: Different processes and corresponding measurements contributing to the determi- 
nation of p, fj, Cb a , 4>B d , Cb 3 ,<Pb s an d C £K . Am K is not considered due to the fact that 
the long distance effects are not well under control. 



in the case of B® — B® mixing we define 



/ 7^0 I t/fulll D0\ 

Cn e 2i ^ B " = - gl cS 1 q1 (a = d s) (1) 

where includes only the SM box diagrams, while H^ 1 includes also the NP contribu- 
tions. 3 In the absence of NP effects, C# = 1 and <pB q — by definition. The experimental 
quantities determined from the B® — B® mixings and listed in Tab. El are related to their 
SM counterparts and the NP parameters by the following relations: 

Amf p = C Bd Am s d M , sin 2/3 exp = sin(2/5 SM + 20 B J , a exp = a SM - <\> Bi , (2) 

in a self-explanatory notation. As far as the K° — K° mixing is concerned, we find 
it convenient to introduce a single parameter which relates the imaginary part of the 
amplitude to the SM one: 

]m[(K°\H^\K°)] 
EK Jm[(K°\H^\K°)] ' 1 ' 

This definition implies in fact a simple relation for the measured value of Sk, 

eT = C eK 4 M (K° - K° mixing) . (4) 

ArriK is not considered because the long distance effects are not well under control. There- 
fore, all NP effects which enter the present analysis are parameterized in terms of three 
real quantities, Cs d , (pB d , and C tK . NP in the B s sector is not considered in this case, due 
to the lack of experimental information, since both Am s and Aqp(B s — ► J/ip(p) are not 
measured yet. 



3.1 New Physics effects in the extraction of a from |AF| = 1 pro- 
cesses 

In principle, the extraction of a from B — > titt, pn, pp decays is affected by NP effects in 
|AF|=1 transitions. Actually, in the presence of NP in the strong b — » d penguins, the 
decay amplitudes for B mesons decaying into tttt, pir and pp are a simple generalization 
of the SM ones (given for example in Eqs. (17) and (18) of Ref. [1]). We assume that 

3 Cb and 4>b parameterize NP effects in the dispersive part of the effective Hamiltonian only. 
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NP modifies significantly only the "penguin" amplitude P without changing its isospin 
quantum numbers (i.e. barring large isospin-breaking NP effects). Then, instead of a 
complex penguin amplitude with vanishing weak phase, we have two independent arbi- 
trary complex penguin amplitudes for B and B decays. For example, the amplitudes of 
B — > 7i7T(pp) can be written as 



A + - 
A + - 
A+ o 

A~° 
A 00 
A 00 



-Te~ ia + Pe 



i5 T[ 



-Te ia + p e -^ p e ldp 
1 

"71 
i 

"71 
1 

"71 
i 

"71 



(T + T c e 



ld Tc 



T c e~ ia e i&T - + P e i4>p e iSp 
T c e ia e i6x ° + p e -^ p e i5p 



(5) 



where T, T c , P and P are real parameters, 5p and <5t c are strong phases, a is the angle of 
the UT, and <pp is an additional weak phase. 

The procedure to extract a is exactly the same as in the SM [1], since we assume 
that NP does not affect the AI = 3/2 amplitudes. However, in the NP fit we lose the 
knowledge of the weak phase of the penguins. In spite of this additional free parameter, 
the experimental information available nowadays is sufficient to constrain a even in the 
presence of NP, as shown in Fig. |21 We take T, T c , P and P to be flatly distributed 
in a range larger than the one determined from the fit, and all the phases to be flatly 
distributed in the [0, 2ir) range. Here and in the following figures, dark (light) areas 
correspond to the 68% (95%) probability region. One should notice that these analyses 
bound a through the quantity it — (3 — 7, where 7 comes from the decay amplitudes 
and j3 from B^ — B^ mixing. Therefore, in the presence of NP effects in the \AF\ = 2 
Hamiltonian, this bound should be regarded as a constraint on ct SM — <pB d (see Eq. |2J). 

The reader might notice a contradiction between the discussion above and the results 
of ref. [15,23], in which it is stated that the NP parameters introduced above can be 
eliminated by a redefinition of the T, T c and P parameters. Explicitly, one has 



A + - 
A+ o 



A 



00 



-T x e 


-ia 


1 




"71 
1 


e~ 


"71 





X 



T cX e- ia + P 



x 



with 



T cX 
Px 



T + 



p e i4>P _ p e ~ 



-Upp 



2i sin a 

p e i4>P _ p e -i<)>P 



(6) 



(7) 



,i<5p 



2% sin a 



p e i(a+cj>p) _ p e -i(a+<j> P ) 



iSp 



2% sin a 
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Figure 2: P.d.f. of a from the combination of is o spin analyses ofirn, pre and pp decay 
modes, including NP effects in the \AF\ = 1 Hamiltonian. 



However, the transformations (J7J) are singular for a. — > 0. This implies that there is no 
limited a-priori range for the X parameters. For this reason, the parameterization in 
eq. (JHJ) is of no use for our purpose. 

3.2 Asl: general considerations and the inclusion of |AF| = 1 
New Physics effects 

One can also add the constraint coming from the CP asymmetry in semileptonic B decays 
Asl, defined as 

= r^ -> e+x) - t(b° -> j-x) 

SL = T(B° -f £+X) + T(B° t-X) ' ( ' 

It has been noted in Ref. [13] that, even though the present experimental bound is not 
precise enough to bound p and fj in the Standard Model, A SL is a crucial ingredient of 
the UT analysis once the formulae are generalized according to Eq. (J2J), since this is the 
only constraint that depends on both Cs d and (pB d - 
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m 


0.1797 ± U.UU17 


n 2 


0.1391 ± 0.0193 


™3 


—0.0012 ± 0.0014 


n 4 


-0.0074 ± 0.0020 


n 5 


0.0020 ±0.0007 


n 6 


1.0116 ±0.0826 


n 7 


0.0455 ±0.0144 


n 8 


-0.0004 ± 0.0046 


n 9 


-0.0714 ±0.0170 


nxo 


-0.0041 ±0.0016 


rin 


-0.3331 ±0.2178 




0.0028 ±0.0101 


ni3 


-0.0036 ±0.0033 











Table 4: Magic numbers for the calculation of A^l. The quoted errors correspond to 
Gaussian distributions. 



where Ti2 and M\ 2 are the absorptive and dispersive parts of the B\ — B\ mixing am- 
plitude. At the leading order, A SL is independent of penguin operators, and therefore it 
is also independent of NP in \AF\ = 1 processes. However, at the NLO, the penguin 
contribution should be taken into account. In the SM, the effect of penguin operators is 
GIM suppressed since their CKM factor is aligned with M 12 : both are proportional to 
{V t lV t d) 2 - This is not true anymore in the presence of NP, so that the effects of penguins 
are amplified beyond the SM and the approximation made in Ref. [13] of neglecting this 
contribution is questionable. For our analysis of Ash, we therefore start from the full NLO 
calculation of Ref. [24], allowing for an additional NP contribution to the penguin term 
in the |AF| = 1 amplitude. This introduces two additional parameters (Cp en and </>p en ), 
encoding NP contributions to the penguin part in analogy to what C Bd and (p Bd do for 
the box contribution. Since the penguin amplitude is 0(a s ) with respect to the leading 
contribution, these parameters introduce a smearing in the theoretical determination of 
Asl. The generalized expression of A^h is given by 

, 2k f . fn± \ ( n 6 B 2 + n 11 \ sin (f3 + 2<fr Bd ) ( . n 7 B 2 + n 12 \ 

SL = -e^r m(2 M" i+ Bi ) % — r + Bl ) 

Cpen / B 2 



sin (J3 + Pe „ + 20 B J -g 5 (n 5 + n l0 ^- ) \ (10) 



where B\ corresponds to the usual B& parameter for B° — B° mixing, B 2 = 0.84 ± 0.07 



(flat) [25], Rf = y (1 — p) 2 + fj 2 is the length of one of the UT sides, k is defined in Ref. [24] 
and the magic numbers Hi are given in Tab. EJ The Standard Model expression can be 
recovered in the limit Cx — > 1 and <px ~ * (where X = Ba, Pen). Eq. (I10J) contains 
NLO QCD and l/m^ corrections; the latter have been estimated using matrix elements 
computed in the vacuum insertion approximation, since lattice results are not available. 

To display the main phenomeno logical consequences of Asl, let us consider a simplified 
formula obtained by setting all magic numbers to their central values, and dropping all 
those smaller than 10~ 2 . In this way we get 

2k f , , / 1.0LB 2 -0.33N sm(3 + 2(b B ,) ( B 2 
A SL ~ — |sin(20 fld )(o.l8 + 1- J VBJ (OAA + 0.05-^ 

+ sin (0 Pen + 2<j) Bd ) C Pen (-0.07)|^| . (11) 
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The SM penguin contribution vanishes at this level of accuracy. It is evident from the 
simplified expression in Eq. (fTTj) that the phase <pB d can induce an order-of-magnitude 
enhancement of Agh relative to the SM, while the penguin phase 0p en can induce correc- 
tions comparable to the SM contribution. To be conservative, for our analysis we varied 
Cp en in the range [0,2] with Pcn G [0, 2ir]. This produces only a minor smearing of the 
dominant effects due to NP in \AB\ = 2 transitions. 

3.3 Results of the analysis and constraints on |AF| = 2 NP con- 
tributions 

To obtain the constraints on NP we extract Cs d , C eK , p and fj with a flat distribution 
in a range much larger than the experimentally allowed region. The phase 4>B d is taken 
to be flatly distributed in the range [0, 7r]. The generated events are weighted using the 
experimental information on | V^b/V^&l, B — > DK decays (7), €k, B — > pp, pn, and mi 
decays (a), B -> J/VKW and B -> D°h° decays [26,27] ((3), and A SL [18], using the 
technique described in ref. [17]. The output p.d.f.'s for Cs d , C €K , Cb a vs. (f)B d , and 7 vs. 
B d are shown in Fig. 01 and the corresponding regions in the p-fj plane are presented in 
Fig.il 

It is important to remark that the constraints coming from the experimental observ- 
ables allow for an increase in the precision on p and fj with respect to the pure tree-level 
determination. This is clear comparing Fig. ^ to Fig. 

To illustrate the impact of each experimental constraint on the analysis, in Fig. 03 
we show the selected regions in the (ps d vs. Cs d and (pB d vs. 7 planes using different 
combinations of constraints. The first row represents the pre-2004 situation, when only 
IV^/V^I) Am a, Ek and sin 2(3 were available, selecting a continuous band for <ps d as a 
function of 7 and a broad region for Cs d - Adding the determination of 7 (second row), 
only four regions in the (f>B d vs. 7 plane survive, two of which overlap in the (f>B d vs. Cs d 
plane. Two of these solutions have values of cos 2(/3 + (j>B d ) and a — (pB d different from the 
SM predictions, and are therefore disfavoured by (cos2/5) exp and by the measurement of 
(2/5) exp from B — > Dh° decays, and by a exp (third and fourth row respectively). On the 
other hand, the third solution has a very large value for Asl and is therefore disfavoured 
by ^4gL P , leading to the final results already presented in Fig. El 

In Tab.Elwe give the numerical results for the NP parameters and some of the relevant 
UT quantities, for which we show the output distributions in Fig. El A comment is 
needed for the case of Am s : the output distribution reported in Fig. [7| represents the SM 
contribution only (i.e. it corresponds to Cb s = 1). Therefore this numerical result should 
not be taken as a prediction for Am s in a general NP scenario in which Cb s 7^ 1- The 
conclusion that we can draw from the output distribution of Am s is most easily read from 
the compatibility plot 4 shown in Fig. [7| a value of Am s > 30 (36) ps _1 would imply the 
presence of NP in B s — B s mixing at the 2 (3) a level. On the other hand, from the similar 
result in the contest of the Standard Model [1] one can still conclude that Am s > 29 (34) 
ps _1 would imply the presence of NP at the 2 (3) a level (but not necessarily in the B s 
sector) . 

Before concluding this section, let us analyze more in detail the results in Fig. |21 
4 The method used to calculate the level of agreement in the compatibility plot is explained in [1]. 
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Figure 3: Output P.d.f.'s for C eK (top-left), Cs d (top- center) , <p Bd (top-right), and 2D 
distributions of 4>B d vs. Cs d (bottom-left) and (j>B d vs. r j . Dark (light) areas correspond to 
the 68% (95% ) probability region. 



Writing 



A o 2i P _l_ A 2i{P+<j) NP ) 

2i</, n . ^SMe M + /iNpe VM V1N± ^ 



^ 6 ' (12) 

and given the p.d.f. for Cs d and we can derive the p.d.f. in the (A^p/Asm) vs. 0np 
plane. The result is reported in Fig. [HI We see that the NP contribution can be substantial 
if its phase is close to the SM phase, while for arbitrary phases its magnitude has to be 
much smaller than the SM one. Notice that, with the latest data, the SM (4>B d — 0) is 
disfavoured at 68% probability due to a slight disagreement between sin 2f3 and |VL&Ac&|- 
This requires Anp ^ and 0np ^ 0. For the same reason, 0np > 90° at 68% probability 
and the plot is not symmetric around 0np = 90°. 

A similar parameterization has been used in ref. [28]. Comparing our Fig. |H] with Fig. 5 
of ref. [28], one notices small differences. However, since they are using the statistical 
method of ref. [15], they are plotting areas corresponding to "at least" confidence levels, 
so that their areas are expected to be larger than ours. 

Assuming that the small but non- vanishing value for (pB d we obtained is just due to a 
statistical fluctuation, the result of our analysis points either towards models in which new 
sources of flavour and CP violation are only present in b — > s transitions, a well-motivated 
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p 

Figure 4: The selected region on p-fj plane obtained from the NP generalized analysis. 
Selected regions corresponding to 68% and 95% probability are shown, together with 95% 
probability regions forj (from DK final states) and \ V u b/V c b\. 



possibility in flavour models and in grand-unified models, or towards models with no new 
source of flavour and CP violation beyond the ones present in the SM (Minimal Flavour 
Violation). This second possibility will be studied in detail in Sectional 

4 Constraints on NP from |AS| = 2 or |AB| = 2 tran- 
sitions only 

A complementary information to the one presented in the previous section is obtained by 
allowing NP contributions to be present only in |A5| = 2 or \AB\ = 2 transitions. This 
can be useful to test models beyond the SM in which NP contributions are expected to 
affect dominantly only one of these two sectors, and is also the starting point to update 
previous analyses of NP in \ AS\ = 2 or \ AB\ = 2 processes in supersymmetry [29,30] or 
in any other given model. 

Allowing NP to affect only C EK , we obtain the results for the UT parameters, for 
C EK and for p and fj reported in Fig. |U] and in Tab. El The determination of the UT is 
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Figure 5: 2D distributions of <pB d vs. Cs d (left) and 4>B d vs. 7 (right) using the following 
constraints: i) \V u b/V c b\, Am^, Ek and sin 2(3 (first row); ii) the constraints in i) plus 
7 (second row); Hi) the constraints in ii) plus cos 2/3 from — > J/ipK* and (3 from 
B — > (third row); iv) the constraints in ii) plus a (fourth row). 
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Figure 6: From iop to bottom and from left to right, the p.d.f. 's for p, fj, a, sin 2(3, 7 and 
2(3 + 7. The red (darker) and the yellow (lighter) zones correspond respectively to 68% 
and 95% of the area. These results are obtained in the presence of NP in all the processes 
entering the UT analysis. 
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Generalized XJTfit analysis in the presence of NP 



Standard Solution (7 > 0) Non-Standard Solution (7 < 0) 



UT parameters 


P 


0.246 ± 0.053 ([0.115,0.370] @95%) 


[-0.230,-0.212] 


@95% 


V 


0.379 ± 0.039 ([0.277,0.463] @95%) 


[-0.398,-0.381] 


@95% 


sin 2(3 


0.799 ± 0.037 ([0.694,0.880] @95%) 


[-0.588,-0.574] 


@95% 


7L°] 


57.1 ±7.8 ([37.9,75.4] @95%) 


[-121.5,-118.4] 


@95% 


«[°] 


96.0 ±8.4 ([78.3,116.5] @95%) 


[-44.5, -40.0] ( 


a95% 


2/3 + 7 [°] 


110.9 ±9.2 ([88.8, 128.6] @95%) 


[-158.5,-153.0] 


095% 


Im X t [xl(T 5 ] 


14.9 ± 1.5 ([11.7,17.6] 


@95%) 




Am s [ps _1 ] 


18.0 ±5.3([8.9, 29.6] € 


m%) 




NP related parameters 


C B d 


1.27 ±0.44 ([0.56,2.51] 


@95%) 






-4.7 ±2.3 ([-9.9,1.0] @95%) 


[39.0,39.8] <§ 


95% 




0.95 ±0.18 ([0.64,1.44] @95%) 


[-0.71,-0.59] <S 


^95% 



Table 5: Results of the NP generalized analysis on UT parameters. The values for CB d , 
4>B d and C €K are reported. The second solution is excluded at 68% probability level so we 
quote the 95% ranges only. 




A mjps' 1 ] A m s [ps' 1 ] 



Figure 7: P.d.f. (left) and compatibility plot (right) for the SM contribution to Am s in 
the presence of NP in all the quantities entering the UT analysis, setting Cb s = 1- 



essentially equivalent to the SM one, since only Ek is missing in this case. 

For the case in which NP only enters Bd — Bd mixing, the results are given in Fig. El 
and in Tab. |H1 The main difference with the results in the previous section is that one 
can use Ek to eliminate the solutions with negative 7. 
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Figure 8: P.d.f. in the (Anp/^sm) vs. 0np plane for NP in the \AB\ = 2 sector (see 
Eq. 



5 Minimal Flavour Violation models 

We now specialize to the case of MFV. Making the basic assumption that the only source 
of flavour and CP violation is in the Yukawa couplings [6], it can be shown that: 

1. The phase of \AB\ = 2 amplitudes is unaffected by NP, and so is the ratio 
Am s / Arrid. This allows the determination of the Universal Unitarity Triangle 
independent on NP effects, based on |V^/Vrf,|, 7, Acp(B — > J/^K^), (3 from 
B -> D°h°, a, and Am s /Am d [5]. 

2. For one-Higgs-doublet models, and for two-Higgs-doublet models at low tan/3, all 
NP effects in the UT analysis amount to a redefinition of the top box contribution 
to \ AF\ = 2 processes S (x t ) — > S (x t ) + SS . 

3. For two-Higgs-doublet models with large tan/3, NP enters in a similar way with re- 
spect to the low tan f3 case, but this time one cannot relate the parameter redefining 
So(x t ) in the B sector to the similar term in the K sector. Therefore, two different 
redefinitions must be made for the B and K sectors: S (x t ) — > S Q (x t ) + SSq' k . 
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Figure 9: The p.d.f. for C £K (left) and the p — fj plane (right). In the left plot, the darker 
and the lighter zones correspond respectively to 68% and 95% of the area. These results 
are obtained in the presence of NP in K° — K° mixing only. 



We perform three different analyses, corresponding to the points 1.-3. above. First, we 
present the determination of the UUT, which is independent of NP contributions (Sec. l5.ljl 
in the context of MFV models. Then we add to the analysis the NP parameter 5 So and 
constrain it, together with p and fj, using also the neutral meson mixing amplitudes. 
Finally, we consider the case SSq ^ SSq and determine the constraints on p, fj and these 
NP parameters. We take 5S , 6~Sq and 5Sq to be flatly distributed in a range much larger 
than the experimentally allowed region. 

5.1 Universal Unitarity Triangle 

In Fig. ^2 we show the allowed region in the p — fj plane for the UUT, and in Fig. ^] 
we plot the p.d.f. 's for several UT quantities. The corresponding values and ranges are 
reported in Tab. The most important differences with respect to the general case are 
that i) the lower bound on Am s forbids the solution in the third quadrant, and ii) the 
constraint from sin 2/? is now effective, so that we are left with a region very similar to 
the SM one (for the reader's convenience, we also report results of the SM UT analysis in 
Tab. EJ). The values in Tab. [7|are the starting point for any study of rare decays and CP 
violation in MFV models. See Ref. [10] for a recent analysis based on the results of this 
work. 

5.2 Constraints on NP contributions in MFV models 

We now determine the allowed ranges of NP contributions to |AF| = 2 processes, both 
in the small and large tan/3 regime. Furthermore, using the conventions of Ref. [6], we 
quantify the scale of NP that can be probed with the UT analysis. 

Let us start by considering MFV models with one Higgs doublet or low/moderate 
tan/3. In this case, all NP effects in \AF\ = 2 transitions are due to the effective Hamil- 
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to 68% and 95% of the area. These results are obtained in the presence of NP in — B% 
mixing only. 
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Generalized XJTftt analysis in the 


presence of NP 




\AS\ = 2 only 


\AB d \ = 2 only 


UT parameters 


P 


0.267 ± 0.056 [0.145,0.368] 


0.246 ± 0.038 [0.166,0.317] 


V 


0.319 ± 0.034 [0.257,0.387] 


0.372 ± 0.028 [0.318,0.424] 


sin 2(3 


0.730 ± 0.028 [0.674,0.781] 


0.794 ± 0.033 [0.727,0.854] 


7[°] 


50 ± 9 [35,69] 


56 ± 5 [46,67] 


«[°] 


107 ±9 [87,123] 


97 ±6 [86,108] 


2/? + 7 [°] 


100 ±10 [80,118] 


109±6[96,120] 


Im A t [xl0~ 5 ] 


12.6 ± 1.4 [10.1,15.4] 


14.9 ± 1.0 [13.0,16.8] 


Am s [ps^ 1 ] 


22.6 ±4.2 [15.2,30.8] 


17.4 ±2.0 [14.8,22.4] 


NP related parameters 




1 


1.25 ± 0.21 [0.84,1.69] 


0B d [°] 





-4.6 ± 2.0 [-8.5,-0.7] 




1.10 ±0.21 [0.73,1.59] 


1 



Table 6: Results of the NP generalized analysis on UT parameters, when only NP contri- 
butions to the \ AS\ = 2 (left) and \ ABd\ = 2 (right) processes are considered. The values 
for Cs d , 4>B d and C 6K are reported. The quoted errors represent 68% [95%] probability 
ranges. 



tonian 

~(Ql\fc1»Ql)\ (13) 

with (\Fc)ij — YtV£V t j for i ^ j and zero otherwise, Y t the top quark Yukawa coupling, 
A the scale of NP and a an unknown (but real) Wilson coefficient. The value of a 
can range from order one for strongly interacting extensions of the SM to much smaller 
values for weakly interacting theories and/or symmetry suppressions analogous to the 
GIM mechanism in the SM. It is now trivial to project this onto the SM |AF| = 2 
effective Hamiltonian: it amounts only to a modification of the top quark contribution 
to box diagrams. Normalizing the NP Wilson coefficient to the SM effective electroweak 
scale 6 A = Y t sin 2 Q w M w ja « 2.4 TeV, we obtain 

S (x t )^S (x t ) + 5S , 5S = 4a^ 2 . (14) 

We can therefore determine simultaneously the shape of the UT and 5 So from the standard 
UT analysis. Then, choosing as reference values a = ±1, we can translate the constraints 
on 5 So into a lower bound on A. We obtain (see Fig. [T3*j) : 

, 3.6 TeV @95% Prob. for positive 5S , s 

A> { (15) 
1 5.1 TeV @95% Prob. for negative 5S 



5 Here and in the rest of this section we follow the notation of Ref. [6]. 

6 i.e. the scale obtained by writing the SM contribution to \ AF\ = 2 transitions in the form of Eq. (|13fl 
with coefficients a of order one. 
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Figure 11: The selected region on p-rj plane obtained from the determination of the UUT. 



Also in this case, we can obtain predictions for UT parameters, together with a con- 
straint on NP contributions (see Tab. IHJ). 

In the case of large tan /3, the situation changes since the bottom Yukawa coupling is 
not negligible anymore, and it can distinguish transitions involving b quarks from those 
involving only light quarks. This spoils the correlation of \AB\ = 2 with |A5| = 2 
amplitudes, so that two uncorrelated parameters 5Sq and SSq are required in this case, 
to take into account NP contributions to Bd, s -Bd, s and K-K mixing. In a global fit, made 
by using all the available inputs, Am^ and Amd/Am s determine the value of SSq, Ek 
fixes SSq , while p and fj are given by the combination of all the other constraints. 

Performing this analysis, we bound the UT parameters as given in Tab. |Hl limiting 
the NP scale to be: 



A > 



A > 



L , from B dtS - B d , s mixing 



2.6 TeV @95% Prob. for positive 8S$ 
4.9 TeV @95% Prob. for negative 8S& 

3.2 TeV @95% Prob. for positive 6S£ r T ^ - 

„ from K — K mixing 
4.9 TeV @95% Prob. for negative 5S£ 



(16) 



The output distributions for SSq and SSq are given in Fig. [T31 
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yf] 2p+7[°] 

Figure 12: From top to bottom and from left to right, the p.d.f. 's for p, fj, a, sin2/? ; 7 and 
2(3 + 7. The red (darker) and the yellow (lighter) zones correspond respectively to 68% 
and 95% of the area. These results are obtained from the UUT analysis. 
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Universal Unitarity Triangle analysis 





UUT (68%) 


UUT (95%) 


SM (68%) 


SM (95%) 


UT parameters 


P 


0.259 ± 0.068 


[0.107, 0.376] 


0.216± 0.036 


[0.143, 0.288] 


fj 


0.320 ± 0.042 


[0.241, 0.399] 


0.342±0.022 


[0.300, 0.385] 


sin 2(3 


0.728 ± 0.031 


[0.668, 0.778] 


0.735± 0.024 


[0.688, 0.781] 


a[°] 


105 ± 11 


[81, 124] 


98.5± 5.7 


[87.1, 109.8] 


J[°] 


51 ± 10 


[33, 75] 


57.6± 5.5 


[46.8, 68.7] 


(2/5 + 7) [°] 


98 ± 12 


[77, 123] 


105.3± 8.1 


[89.6, 121.4] 


Im A t [xl0~ 5 ] 


12.7 ± 1.7 


[9.7, 15.9] 


13.5 ± 0.8 


[12.0, 15.0] 


Am s [ps -1 ] 


20.6 ± 5.6 


[10.6, 32.6] 


20.0 ± 1.8 


[15.5, 24.2] 



Table 7: Results of the UUT analysis. For convenience, the SM results from Ref. [1] are 
also reported. 



Minimal Flavour Violation analysis 





low/moderate tan/3 


large 


tan (3 




68% 


95% 


68% 


95% 


p 


0.216 ± 0.058 


[0.109, 0.361] 


0.231± 0.067 


[0.112, 0.375] 


V 


0.351 ± 0.032 


[0.265, 0.406] 


0.347±0.036 


[0.254, 0.404] 


sin 2(3 


0.733 ± 0.027 


[0.679, 0.781] 


0.731± 0.027 


[0.673, 0.781] 




98.6 ± 9.5 


[81.6, 121.7] 


101± 11 


[82, 124] 


7[°] 


57.6 ± 9.1 


[35.7, 79.1] 


55± 11 


[34, 74] 


W + 7) [°] 


104 ± 10 


[80, 122] 


102± 12 


[77, 121] 


Im A t [xlO -5 ] 


13.6 ± 1.4 


[10.1, 16.0] 


13.4 ± 1.9 


[9.7, 16.3] 


Am s [ps -1 ] 


19.5 ± 2.6 


[15.0, 31.7] 


22.6 ± 5.4 


[15.5, 35.1] 



Table 8: Results for UT parameters from the MFV generalized analysis. 



It is instructive to observe the two-dimensional plot of 5S$ vs. 5Sq in Fig. ED within 
models with only one Higgs doublet or with small tan (3, the two <5's are bound to lie on 
the line 5S$ = 5Sq. The correlation coefficient R provides a measure of this relation. 
We find R = 0.52 giving no compelling indication on the value of tan/3. 

6 Model Independent constraints on New Physics in 
the |AF|=2 sector in year 2010 

We present an exercise on the knowledge of the UT parameters within the generalized 
NP analysis in a possible scenario in year 2010. At this date, the B factories will have 
completed their data analysis and the LHCb experiment will have started running. 
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Figure 13: P.gL/. of 5 So (top-left), 5Sq vs 5Sq (top-right), 5Sq (bottom-left) and 5Sq 
(bottom-right). See the text for details. 



For this exercise we have assumed a total integrated luminosity for the B factories of 
2 ab _1 and two years of data taking at LHCb, with an integrated luminosity of 4 fb _1 . 
At that time the lattice community will have produced the final numbers from the Tera- 
Flops machines. The 2010 projected values and errors for the quantities which are most 
relevant in UT analysis are given in Tab.|Hl Lattice parameters are taken from [31], while 
the extrapolation of the errors on the experimental measurements is taken from [32-35]. In 
addition to the improvements of existing measurements, we have added new measurements 
in the B s sector from LHCb. In particular the determination of 



Amf p 
sin 2Xs Xp 
( 7 - 2 X ,) cxp 



C B Am^ M 
sin(2 X ^ M + 2<p Bs 



7 



SM 



2X 



SM 



from B s 
from B s 
from 5, 



D s 7l , 

D X K . 



The central values for the different observables have been generated in the SM starting 
from an arbitrarily chosen value of p and fj, so that they are all compatible with each 
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other and the result of the fit is fully "SM like". The reason for this procedure is to 
investigate whether, in the "worst case" scenario of perfect confirmation of the SM, one can 
asymptotically reduce the errors on the NP related quantities introduced in the previous 
sections, and translate the derived constraint into a lower limit on the energy scale for 
NP particles. We start from a picture of what the UT analysis should look like in 2010. 
In Fig. El we show the selected region in the p-fj plane, while in the third column of 
Tab. El we quote the uncertainties on the various quantities from the UT analysis in 
the Standard Model. This should be taken as a reference for the approaches beyond the 
Standard Model that follow. 

Moving from the SM analysis to the model independent approach of Sec. |HJ we expect 
in the future a sizable improvement of the knowledge of the NP parameters: 

C Bd = 0.98 ±0.14 (j) Bd = (-0.1 ± 1.3)° 

C Ba = 0.99 ±0.12 <f> Ba = (0.0 ± 1.3)° 

C eK = 1.00 ± 0.10 (17) 

as shown in Fig. E3 

In the same future scenario, one can repeat the MFV analysis, both determining UT 
parameters using the UUT approach and adding the information from (NP sensitive) mix- 
ing quantities to bound the NP scale. The expected errors on the relevant UT quantities 
are summarized in the forth and fifth columns of Tab. El 

If no evidence of violation of the Standard Model will emerge from B physics in the era 
of direct NP search at LHC, this generalized approach will replace the present UT analysis 
as the default procedure. So, it is important to remark the fact that the generalization of 
the analysis costs an increasing of about 10% of the errors, which is not a huge price to 
pay if compared to the gain in terms of the larger physics scenario. 7 In Fig. E3 we give 
a hint of what the UUT analysis would look like in 2010. 

7 Of course, the output error on Am s is also affected by the absence of Am^ in the fit. 
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Observable 


projected value & error 


sin OR 

bill ZifJ 


u.uyo nz u.uiu ^i.ft/oj 


n/\°] 


1 D4 4- ^ 


7L J \ ury ) 


54 4- 5 




n QQn 4- n nzL7 


i , B s \jB Bs [MeVJ 


U.Z(0 ± U.U14 [p/Q) 


r 


l.zUU ± U.Uof {o/oj 


\Vcfj -(mcl+excl) (10 J 


41.7 ± (J. 4 (0.97o) 


|l4 fe |-(mci+excl) (10 *) 


do. 4 ± l.o (4.27oJ 


Am d [ps X J 


0.503 ± 0.003 (0.6%) 


,™ r .\ n 
m t [IjeVJ 


1 11 ± O.U 


A 


0.2240 ± 0.0008 


Am s [ps -1 ] 


20.5 ± 0.3 


sin2x s (J/^0) 


0.031 ± 0.045 


(7-2 X ,)[°] (D S K) 


51 ± 10 



Table 9: Projected values and errors in year 2010 for the most relevant quantities entering 
in the UT analysis. The central values are chosen such that the constraints are perfectly 
compatible within the SM. 



In this framework, one should expect to increase the lower bound on A when the NP 
sensitive quantities are added to the UUT fit. To have a quantitative example of the 
expected improvement, we used the input listed above for the 2010 scenario, obtaining 
the SSo distributions shown in Fig.[T7| From these distributions, we get A > 7.5(6.6) TeV 
at 95% probability, in the case of positive (negative) value of 5So, in the case of MFV 
models with one Higgs doublet or low/moderate tan/5. For the case of large tan/?, we get 



A > 



A > 



6.0 TeV @95% Prob. for positive 6S$ 
6.6 TeV @95% Prob. for negative 5S^ 



from B dyS - B, 



i djS mixing 



6.8 TeV @95% Prob. for positive 5S? r T ^ - 

„ from K — K mixing 
5.1 TeV @95% Prob. for negative 5S£ 



18) 



7 Conclusions 



We have performed a model-independent analysis of the UT in general extensions of 
the SM with loop-mediated contributions to FCNC processes. Going beyond the pure 
tree-level determination of the UT already presented in Ref. [1], we have shown how the 
recent measurements performed at B factories allow for a simultaneous determination of 
the CKM parameters together with the NP contributions to |AF| = 2 processes. We have 
found strong constraints on NP contributions that can be as large as the SM ones only if 
the SM and NP amplitudes have the same weak phase. 
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Figure 15: Projected situation in year 2010. p.d.f. distributions for NP parameters: from 
top to bottom and from left to right, distributions of Cs d and <pB d , distributions of Cb s 
and 4>b 3 ! 2D distributions of <pB d vs. Cg. and <pB a vs. Cb s , and distribution of C eK . Dark 
(light) area corresponds to the 68% (95%) probability region. 



Motivated by this result, which points towards models with MFV, we have analyzed 
in detail the UUT. By putting together all the available information, it is possible to 
determine the UT parameters almost as accurately as in the SM case and to constrain the 
additional NP parameters. In this way, we probe dimension-six operators up to a scale of 
5 TeV, to be compared with the SM reference scale of 2.4 TeV and to the sensitivity of 
other rare processes, which reaches scales of 9-12 TeV in the case of b —>■ sj [6]. 

Finally, we have presented a possible scenario for the UT analysis in five years from 
now, taking into account foreseeable progress in theory and experiment, under the pes- 
simistic assumption that the SM perfectly agrees with the data. This exercise allows us 
to assess the sensitivity to NP that we can expect in the near future. The impressive 
accuracy we can reach in this kind of analyses shows the great potential of flavour studies 
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UT analysis in 


zUlU 




Observable 


Input error 


Output 


error 






SM UT 


T TT "TTT"! 

UUT 


P 




0.015 


0.021 


fj 




U.UU ( 


U.U1U 


sin 2/3 


0.010 


0.009 


0.009 


«[°] 


5 


2.1 


3.5 


7[°] 


5 


2.0 


3.2 


2/3 + t[°] 




2.3 


3.4 



Table 10: Projected uncertainties in year 2010 for the SM UT and the UUT analysis, 
using all the inputs of Tab. 
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Figure 16: The selected region on p-fj plane obtained from the UUT analysis in the "year 
2010" scenario. 



in investigating the structure of NP. 
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Figure 17: P.d.f. of 5 So from MFV fit in the 2010 scenario, in the case of small (left) 
and large values of tan (3, from B^ s -B^ 8 (center) and K-K (right) mixing. 
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Note Added 

While completing this work, we became 
similar to the one performed in Sec. El 



aware of ref. [28], which contains an analysis 
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